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We provide extended evidence that mode-coupling theory (MCT) of supercooled liquids for the 
F12 schematic model admits a microscopic realization based on facilitated spin models with tunable 
facilitation. Depending on the facilitation strength, one observes two distinct dynamical glass tran- 
sition lines-continuous and discontinuous-merging at a dynamical tricritical-like point with critical 
decay exponents consistently related by MCT predictions. The mechanisms of dynamical arrest can 
be naturally interpreted in geometrical terms: the discontinuous and continuous transitions corre- 
spond to bootstrap and standard percolation processes, in which the incipient spanning cluster of 
frozen spins forms either a compact or a fractal structure, respectively. Our cooperative dynamical 
facilitation picture of glassy behavior is complementary to the one based on disordered systems 
and can account for higher-order singularity scenarios in the absence of a finite temperature ther- 
modynamic glass transition. We briefly comment on the relevance of our results to finite spatial 
dimensions and to the F13 schematic model. 
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I. INTRODUCTION 

The glassy state of matter remains an active area of 
research despite decades of study. From a theoretical 
viewpoint, the most fundamental issue concerns the very 
nature of the vitrification process, that is whether the 
amorphous state represents a genuine thermodynamic 
phase or rather a purely dynamical accident. It is well 
know that the issue is made especially difficult by the 
lack of a suitable observable allowing the unambiguous 
identification of the elusive amorphous order, and the 
exceedingly long equilibration times involved in exper- 
iments and simulations. The situation is further com- 
plicated by the fact that theoretical modelling of glassy 
systems has made clear that slow relaxation phenomena 
are ubiquitous and may result from very distinct micro- 
scopic mechanisms [l|. Thus, it may appear rather un- 
likely that a single theoretical framework encompasses 
the large variety of behaviors observed in glassy materi- 
als. For this reason, one is often focused on the attempt 
to predict those features that are thought to be universal 
in a minimal setting, which is the basic philosophy we 
adopt here. 

Mode-coupling theory (MCT), is considered by many 
the most comprehensive first-principle approach to the 
dynamics of supercooled liquids (for reviews, see 043)- 
Despite some limitations, it has been much successful in 
predicting peculiar higher-order singularities which have 
been confirmed by experiments and computer simula- 
tions 

[1-Q3. A key tenet of MCT is that the glass tran- 
sition is a purely dynamical phenomenon unrelated to 
any thermodynamic singularity [TB4I3- This prediction 
is rather puzzling from a statistical mechanics point of 



view because diverging relaxation time-scales are typi- 
cally associated with diverging lengthscales [l8[ , and the 
latter are a consequence of a thermodynamical phase 
transition, as clearly exemplified by traditional critical 
phenomena. It has been argued that this unusual sit- 
uation stems from uncontrolled approximations, which 
crucially miss thermally activated hopping processes re- 
sponsible for low-temperature equilibration. Therefore, 
the dynamical nature of the glass transition sin gula rit y 
has been much debated since its first appearance [lJJ, [2(| 
and much work has been devoted to clarify the status of 
MCT UMl. 



Fruitful developments in the dynamics of mean-field 
disordered systems have subsequently shown the inti- 
mate connection of some spin-glass models with MCT 
and led to the formulation of the so-called random first- 
order transition (RFOT) approach pHl-E^. The most 
important ingredient brought about into the discussion 
is the existence of an extra thermodynamic glass phase 
at a Kauzmann temperature T K below the dynamical 
glass transition predicted by MCT. This has immediately 
suggested that while the dynamical transition is an arti- 
fact of mean-field approximation (and therefore doomed 
to disappear in real systems), genuine glassy behavior 
would be a more profound reflection of an underlying 
non-trivial Gibbs measure corresponding to a one-step 
replica-symmetry breaking spin-glass phase. Such de- 
velopments have been generalized in several ways and 
have led to a low-temperature extension of MCT to off- 
equilibrium situations, revealing new interesting features 
such as aging phenomena [32h35| and effective tempera- 
ture [36l. f37(| . 

The thermodynamic RFOT perspective is in appar- 
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ent conflict with a different line of thought inspired by 
facilitated (or kinetically constrained) models first intro- 
duced by Fredrickson and Andersen [38[ (see Ref. [39| for 
a review). These models have a trivial thermodynam- 
ics, and so their glassy behavior has a purely dynamical 
origin by construction. For this reason, they appear to 
be at odds with any thermodynamic perspective on the 
glass transition problem. Moreover, since the simplest 
(so-called non-cooperative) version of such systems does 
not have neither a thermodynamic nor a dynamical tran- 
sition, the facilitation approach also appears to be incon- 
sistent with MCT. Therefore, the notion of dynamical 
facilitation seems to provide an alternative way to the 
glass formation quite in contrast to the MCT and RFOT 
approach (a viewpoint advocated, e.g., in Ref. [40(). It 
has long been noticed, however, that several peculiar fea- 
tures of MCT and RFOT are also present in facilitated 
and kinetically constrained systems (4l| - |43| . In fact, we 
have recently shown that one can construct a class of fa- 
cilitated spin systems that actually reproduces the two 
most relevant MCT scenarios 0, EE} . This means that 
the dynamic facilitation approach on one hand, and the 
MCT and RFOT scenario on the other, should not be 
considered as alternative but rather as complementary 
representations of glassy dynamics. Our statement is 
consistent with recent results showing that, in a specific 
case, there is an exact mapping between facilitated and 
disordered spin systems pa ]. 

In this paper we present detailed numerical results sup- 
porting the idea that cooperative facilitated systems pro- 
vide a microscopic realization of MCT including higher- 
order bifurcation singularity. The rest of the paper is 
organized as follows. To set the stage, in the section 2, 
we briefly review the basic predictions of MCT approach 
to glass transition. In the section 3, we introduce facil- 
itated spin systems and recall their connections to the 
bootstrap percolation problem on a Bethe lattice. In the 
section 4 we define the model we study in this paper and 
report the exact results for the phase diagram and the 
arrested part of correlations. Numerical simulation for 
the relaxation time and the critical decay law exponents 
are presented and compared with MCT predictions in the 
section 5. Finally, in the section 6 we present the conclu- 
sions and possible directions of future works including a 
brief discussion on the possible relevance of our results to 
finite spatial dimension and to the F13 schematic model. 
Some technical calculations are reported in the appendix. 



II. MODE-COUPLING THEORY: THE F i2 
SCENARIO 

A comprehensive description of MCT can be found 
in Refs. Q. Here, we outline only those aspects which 
are pertinent to the present work and that are useful 
for comparison with our theoretical and simulation re- 
sults. In particular, since we focus on schematic models 
we shall disregard the wavevector dependence of the rel- 



evant physical quantities. In the MCT, the microscopic 
system dynamics is projected, through the Zwanzig-Mori 
formalism, on a set of observables which is thought to 
be physically relevant, namely the local density fluctu- 
ations, Sp(t). The correlator of these quantity (p(t) = 
(5p(t)Sp(0)) / (\Sp\ 2 ) is shown to exactly obey the evolu- 
tion equation: 

<j){t) + T 4>{i) + f m(t-s) <]>(s)ds = 0, (1) 
Jo 

where To is the characteristic microscopic timescale of 
the system and we assume, for simplicity, overdamped 
local motion appropriate e.g. to hard-sphere colloids (for 
molecular liquids one should add an inertial term). The 
memory kernel m(t — s) takes into account the retarded 
friction effect which arises from the caging of a particle 
by its neighbors. It is controlled by the static structure 
factor S(q) = (l/N)(\Sp\ 2 ) and does not involve any ther- 
modynamic singularity. In the dilute limit the memory 
function is negligible and the relaxation is exponential. 
At high-density or low-temperature, the memory kernel 
cause a viscosity increase through a feedback mechanism 
that eventually leads to the structural arrest of fluctu- 
ations. For a general system, the explicit expression of 
m{t — s) is rather complicated, and in order to obtain 
quantitative predictions, one is forced to introduce some 
approximations. In schematic models the memory func- 
tion is approximated by a low order polynomial of cor- 
relators, with coupling constants that depend solely on 
S(q). A key quantity in describing the liquid-glass transi- 
tion is the so-called nonergodicity parameter (also known 
as Edwards- Anderson parameter in the spin-glass litera- 
ture), that is the long-time limit of the correlator, 

$ = lim <j)(t). (2) 

t— >oo 

In the fluid state, the system is ergodic and 4>(t) decays 
to zero with time, $ = 0. When $ is finite the system 
is unable to fully relax, its dynamical is arrested and the 
liquid becomes a glass. For schematic models one can 
easily derive the liquid-glass phase diagram. Taking the 
long-time limit of Eq. (JTJ, one obtains the bifurcation 
equation: 



whose solutions possibly give different types of ergodic- 
nonergodic transitions. Of interest here is the so-called 
F12 schematic model which is defined by the memory 
function 

m[4>(t)] = v 1 <j)(t)+v 2 (j)' 2 (t), (4) 

where v± , vi are the coupling parameters controlling the 
system state. In this case, by solving the bifurcation 
equation one finds a discontinuous liquid-glass transition 
line, Vi = \/4:V2 — i>2 with vi G [1,4], at which $ jumps 
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from zero to a finite value $ c , with a square- root singu- 
larity, $ — $ c ~ e 1 / 2 , where e is the distance from the 
transition line. For V2 € [0, 1] one finds a continuous 
liquid-glass transition, vf = 1, across which $ smoothly 
departs from zero in a linear fashion, $ ~ e. 

On the liquid side of the relaxation dynamics, MCT 
makes several specific predictions. Near the discontin- 
uous transition line, the relaxation has a peculiar two- 
step form. The approach and departure from the critical 
plateau at $ c , respectively defining the /? and a relax- 
ation regimes, are described by the so-called critical de- 
cay laws: 



it/re)-", 



p-regime; 
a-regime. 



(5) 



The two characteristic times r and r fl , respectively asso- 
ciated with the a and (3 regimes, increase as power-laws 
near the critical line: 

r ~ t~\ t, ~ e- 1 / 2 ". (6) 

The exponents a, b and 7 are not independent but obey 
the relations: 



A = 



r 2 (i- a ) _ r 2 (i + &) 
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r(l-2a) r(l + 26) 

where T is the Euler's gamma function. The cor- 
relator presents universal scaling (the so-called time- 
temperature superposition principle) in the late a- 
regime: rescaling time by the structural relaxation time 
t, one obtains that data for correlator near the critical 
line can be collapsed onto a single master function. 

Near the continuous transition line there is no a-regime 
and one finds a single step relaxation with a weak long 
time tail described by the critical decay law of the /3- 
regime, i.e. the first of Eq. ([5]) with $ c = remains valid. 
However, the exponent characterizing the relaxation time 
is twice that of the discontinuous transition: 

t> ~ e" 1 /". (8) 

It is important to notice that the exponent a (and there- 
fore b and 7) depends on the actual location of the criti- 
cal point along the glass transition lines. So this provide 
us with the possibility of carefully checking the peculiar 
MCT prediction connecting the exponents of the criti- 
cal decay law and the relaxation time near the transition 
over a wide range of parameters. In particular, when 
the point (u*,^) = (1,1) separating the discontinuous 
and continuous transition is approached, a tends to van- 
ish. Near this glass singularity, which is also called a 
degenerate A3 point, the power-law behavior turns into 
logarithmic relaxation. 



III. FACILITATED SPIN SYSTEMS AND 
BOOTSTRAP PERCOLATION ON BETHE 
LATTICE 

The basic assumption of the dynamical facilitation ap- 
proach is that, on a suitable coarse-grained lengthscale, 



one can model the structure of a liquid by an assembly 
of higher/lower density mesoscopic cells with no ener- 
getic interaction. A binary spin variable is assigned to 
every cell depending on its solid- or liquid-like structure. 
The next crucial step is to postulate that there exists 
a timescale over which the effective microscopic dynam- 
ics can be mapped onto a 'simple' form: local changes 
in cells structure occur if and only if there is a suffi- 
ciently large number, say /, of nearby liquid-like cells (/ 
is called the facilitation parameter). The latter assump- 
tion is admittedly quite remote from the actual liquid 
dynamics, and indeed very difficult to derive by analyt- 
ical means. Nevertheless, it can be justified on phys- 
ical grounds: it mimics the cage effect and gives rise 
to a large variety of remarkable, and sometimes unex- 
pected, glassy features, even if the thermodynamics is 
completely trivial. They include stretched exponential 
relaxation, super- Arrhenius equilibration time, physical 
aging, effective temperature, exponential multitude of 
blocked states, dynamical heterogeneity, etc. 

Facilitated spin models consist of N non-interacting 
spins Cj = ±1, i — 1, . . . , N with Hamiltonian 



N 



(9) 



Spins evolve according to a Metropolis-like dynamics: at 
each time step a randomly chosen spin is flipped with 
transition probability: 



W <7„; -> 



-Cj) 



iin{l,e-' UT ' /feBT }, (10) 



if and only if at least / of its z neighboring spins are in 
the state —1. When the temperature T is low enough, 
the fraction of —1 spins is exponentially small, e ~ h / k sT ^ 
therefore spin flipping is rare and sluggish relaxation en- 
sues. We consider a Bethe lattice with fixed connectivity, 
specifically, a regular fc-graph with coordination number 
z = k + 1: besides one branch going down, each node has 
further k branches going up. The tree-like, local struc- 
ture of the Bethe lattice allows for some exact calcula- 
tions concerning the asymptotic dynamics. Let us call 
B the probability that, without taking advantage of the 
configuration on the bottom, the spin o~i is in, or can 
be brought to, the state -1 by only rearranging the sites 
above it. B verifies the fixed-point equation: 

B = {l- P )+ P Y j ( k \B k -\l-BY (11) 

where p is the probability that the spin is +1 in thermal 
equilibrium, i.e.: 



1 



P 



l + e-VT' 



(12) 



where we set h/k B = 1 hereafter. Eq. ([TTj) is cl osely 
related to that of bootstrap percolation (BP) [43|, l47j . 



4 



a problem that is known to emerge in a wide vari- 
ety of contexts (for a review, see Refs. [HI, H^|). In 
BP every lattice site is first occupied by a particle 
at random with probability p. Then, one randomly 
removes particles which have less than m neighbors. 
Iterating this procedure leads to two possible asymptotic 
results 47]. If the initial particle density is larger than 
a threshold p c there is an infinite cluster of particles 
that survives the culling procedure, whereas for p <Pc 
the density of residual particles is zero. In Ref. (47[ 
it has been shown that the probability 1 — B that a 
particle is blocked because it has at least m neighboring 
particles above it which are blocked (without taking 
advantage of vacancies below) satisfies the self-consistent 
Eq. (fTT|) . By exploiting the relation with the BP one 
can characterize the dynamics on a Bethe lattice. There 
are three qualitatively distinct cases that we now discuss. 



A. Noncooperative dynamics 

For / = 1 , any spin only needs one nearby down-spin to 
flip and the dynamics is called noncooperative. Eq. 
becomes 



1 - B 



(13) 



obviously implying that there is only the solution B = 1 
for any p, i.e. with probability one each spin can be 
brought into the up-state through a finite number of al- 
lowed spin- flips. One can see that relaxation dynamics 
has a single-step form with an Arrhenius equilibration 
time at any temperature. Thus, the noncooperative case 
may be useful for modelling strong glasses. For this rea- 
son it will not be considered here. For more complex 
forms of glassy dynamics we need local relaxation events 
involving two or more facilitating (down) spins, i.e. a 
cooperative dynamics. 



B. Cooperative dynamics 

1. Continuous transition 

For / = k or k + 1 there is an additional solution, 
with B < 1, of Eq. (jTTJ) which is obtained when p is 
large enough. The transition to this regime is continuous, 
because for arbitrary small 1 — B, Eq. (jTTJ) becomes 



1- B = p 



f-l 



(1-B 



(14) 



which implies p c — 1/k for cither / = k or k + 1. No- 
tice that in the later case, the equation has the exact 
same form as for / = k. In fact, these cases are com- 
pletely equivalent to conventional percolation ^7|, and 
for this reason, near p c the incipient cluster of perma- 
nently frozen spins (i.e. spins unable to flip because sur- 
rounded by more than / neighboring up spins) has a 



fractal structure. The facilitated spin dynamics of these 
systems has never been explored to our knowledge. We 
notice that the critical temperature in this case is nega- 
tive 1/T C = — ln(fc— 1). This means that in order to study 
relaxation dynamics one must start with initial configu- 
rations having a majority of down spins. 



2. Discontinuous transition 

For 1 < / < k the power of 1 — B on the right and left 
hand side of Eq. (fT4]) are different and therefore only a 
discontinuous transition -at a critical value p c that de- 
pends on k and /- is possible. This transition is charac- 
terized by a square root singularity coming from low tem- 
peratures and is rather peculiar as it produces an unusual 
mixed or hybrid behavior: the fraction of frozen spins 
jumps from zero to a finite value, with divergent fluc- 
tuations as in critical phenomena. This is much similar 
to the behavior of the non-ergodicity parameter in MCT 
and this form of dynamical arrest corresponds to the sud- 
den emergence of a giant cluster of frozen spins with com- 
pact structure. The geometric origin of this behavior has 
been understood quite in detail as being related to the 
divergence of the size of corona clusters near the tran- 
sition [HO, [HI- Several results (42j, including those re- 
lated to large scale cooperative rearrangements responsi- 
ble for slow dynamics [52| , and the TAP-like organization 
of blocked states below the threshold [43( , have already 
suggested a close analogy with MCT and mean-field dis- 
ordered systems. 



IV. FACILITATED SPIN MIXTURES WITH 
TUNABLE FACILITATION: EXACT RESULTS 

The systems considered so far have a fixed facilitation 
determined by the integer number /, and so they can 
exhibit either a discontinuous or a continuous glass tran- 
sition. In this section we consider a recently introduced 
class of facilitated spin systems in which the facilitation 
strength can be smoothly tuned [44j . This is obtained by 
making the facilitation parameter a lattice site dependent 
random variable, in close analogy with the random boot- 
strap percolation problem studied by Branco [53| (see 
also Refs. [H], HH for recent developments in the con- 
text of complex networks). This means that for every 
site i the facilitation of the spin <jj is a random variable 
generally described by the probability distribution 

fe+l 

where the weights {wi} controlling the facilitation 
strength satisfy the condition < wg < 1 and the nor- 
malization 



fc+i 
e=o 



1. 



(16) 
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By suitably changing the coefficients {wi} one can thus 
explore a variety of different situations, e.g. the robust- 
ness of the glass phase against dilution and, more in- 
terestingly, the crossover between the discontinuous and 
continuous glass transition or the existence of multiple 
glassy states. We remark that mixtures of particles with 
different kinetic constraints were first considered in a 
granular matter context, in the attempt to understand 
particle segregation phenomena under gravity by a purely 
dynamical approach (56l - [58j ). Hereafter, we shall denote 
with {■ ■ ■) f the average over the probability distribution 
Eq. ([151). Following Refs. H3, Ezl H3 , we can easily gen- 



eralize the results presented in the previous section and 
evaluate the probability B as: 

B = l-p + p|EP^"(l^)j ■ (17) 

As we have already anticipated, the precise nature of the 
glass transition depends on the behavior of the fraction 
of permanently frozen spins that, in this framework, can 
be exactly computed from B as follows 



(18) 



where 



If -2 

£ 

\m=0 



B m (l-B) 



(19) 



The two contributions in Eq. (1181) represent the proba- 
bility that a spin is frozen in the ±1 state, respectively. 
In order to evaluate $, in fact, one must take into ac- 
count that in the neighborhood of both up and down spin 
states, there must be less then / facilitated spins. Let us 
consider, for clarity, the homogeneous facilitated case (for 
heterogeneous facilitation, one just needs to further av- 
erage over the distribution Eq. (|T5j) b To be blocked, an 
up spin (found with probability p) may have, among its 
z neighbors, up to / — 1 down spins. These neighbors 
already have one spin up (the main spin) and the prob- 
ability that each one is facilitated is precisely B. The 
overall contribution of all permutations is thus given by 
the first term in the above equation. When the main spin 
is down (with probability 1 — p), we write a term similar 
to the first one, but replacing B with 1 — h. h is the 
probability of a neighbor of the main spin (that is down) 
being up (with probability p) and blocked (because be- 
yond the main spin, we consider the possibility of only 
/ — 2 further down spins). Notice that now to evaluate 
h we consider the possible permutations of / — 2 down 
spins among the 2 — 1 neighbors. When the facilitation 
is no longer homogeneous, we should average both $ and 
h over the distribution V(fi). 

The calculation of the phase diagram and the fraction 
of frozen spins for different types of mixtures on a Bethe 
lattice with z — 4 is detailed in the appendix. Here, 
we shall focus on a ternary mixture with facilitations fi 
chosen from the probability distribution 

V(fi) = (!-«) Sf u k-i + (q-r) 8 h . k + r 8 fi , k+1 , (20) 
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FIG. 1: Phase diagram for jz = 4 and facilitation as in 
Eq. (|20[) . The dark region is the glassy phase. The dot- 
ted/solid line is the hybrid/continuous transition. Inset: 
Fraction of frozen spins, Eq. (|A8|) . vs temperature for sev- 
eral values of q and r = 10 -3 . Below q = 1/2, $ jumps to 
a finite value at the transition which is represented by the 
dotted line. 



with < r < q < 1. The phase diagram for a Bethe 
lattice with z — 4 is depicted in Fig. [1] and comprises 
two glass transition lines, of discontinuous and continu- 
ous nature, given respectively by 



Tc(q) 



In I 



hi- 



12q) , < q < 



3q-l 



1 2 

2^<3' 



(21) 
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The two lines join smoothly at q = 1/2, corresponding to 
the dynamical tricritical temperature 

T tric = l/ln2. (22) 

The critical temperature depends only on the total frac- 
tion q of / = z—1 and z overconstrained spins, but not on 
their relative importance r. As expected, the glass tran- 
sition temperature increases with q since the constraints 
become stronger and the spins more easily blocked when 
compared to the pure case q = 0. Notice that for q = 2/3 
the critical temperature diverges signalling that there al- 
ways exists a finite fraction of permanently frozen spins 
at any positive temperature. Of course, the above phase 
diagram can be formally extended to negative tempera- 
ture, i.e. for 2/3 < q < 1, by contemplating the possi- 
bility of having an equilibrium distribution of —1 spins 
larger than 1/2. The highest temperature is then — oo 
and corresponds to a fully ordered spins configuration in 
the -1 state. This is the initial configuration one has to 
start with to properly perform simulations of equilibrium 
dynamics at negative temperature. Such a situation will 
not be considered here. The precise nature of the two 
glass transitions in Eq. (|2"Tj) depends on the behavior of 
the fraction of permanently frozen spins $ near T c which, 
in the specific case z = 4, is given in the appendix IA II 
The behavior of $ as a function of the temperature is 
shown in the inset of Fig. [TJ for several values of q and 
with r = 10~ 3 . The height of the critical plateau, shown 
as a dotted line in the figure, decreases as q increases, 
going to as q — > 1/2. We remark that the fraction of 
frozen spins generally depends on r while the phase di- 
agram does not. For 0<g<l/2,$ jumps to a finite 
value $ c = &(T C ) on the transition line, meaning that 
the infinite cluster of frozen spins has a compact struc- 
ture. This structure turns out to be quite resilient against 
random inhomogeneities in the facilitation strength, in a 
rather large range of value of q. The critical exponent 
associated to the order parameter $ is obtained by 
expanding the above equations in the small parameter 
e = T — T c . As expected, when $ c has a finite jump, 
we find $ — <I>c ~ ^ i with = 1/2, the typical square- 
root dependence well known in MCT and in other sys- 
tems with a hybrid transition. For 1/2 < q < 2/3 the 
glass transition changes nature as $ departs smoothly 
from zero with a power-law behavior $ ~ and the cor- 
respondence with standard percolation implies that the 
giant cluster of frozen spins has a fractal structure. In- 
terestingly, in this case (3 depends on r: for r — one has 
= 2, while as soon as a negligibly small amount of spins 
with fi = 4 is introduced into the system, r > 0, one has 
/3 = 1. The robustness of the latter behavior reflects the 
fact that the mass of the fractal cluster of frozen spins is 
essentially dominated by the dangling ends, that is, those 
parts of the cluster which are connected to the backbone 
by a single frozen spin. Only when r = 0, dangling ends 
are completely removed from the infinite cluster of frozen 
spins, and the exponent changes to = 2 [53j]. Hence, the 
general scenario emerging for an arbitrary ternary mix- 



ture is that of two distinct glass transitions with 0=1/2 
(for the discontinuous case), and = 1 (for the continu- 
ous one). These critical exponents reproduce exactly the 
MCT results for the F12 schematic model Q. Clearly, 
it is important to establish how these exponents are at- 
tained when the critical line is approached, Also, it is 
by no means obvious that systems sharing the same or- 
der parameter exponents in the glassy phase will also 
have identical critical dynamical exponents. This type of 
problems will be addressed in the next section. 



V. NUMERICAL SIMULATIONS OF A SPIN 
MIXTURE WITH DISCONTINUOUS AND 
CONTINUOUS GLASS TRANSITIONS 

Since the models we have introduced in the previous 
section are only partially solvable, we now turn to nu- 
merical simulation. This will give us the possibility of 
verifying the power-law behavior of relaxation time with 
temperature and of critical decay with time and, most im- 
portantly, the peculiar prediction of MCT that connects 
the exponents of these two laws near the continuous and 
discontinuous transition line. One obvious advantage of 
working with facilitated spin models is that one can im- 
mediately overcome the problem of the long simulation 
times needed for thermal equilibration. Moreover, the ex- 
act knowledge of both the critical temperature and the 
arrested part of correlations makes the procedure for es- 
timating power-law exponents much more reliable, as it 
decreases the number of fitting parameters from three to 
one. In fact, in real systems, for a given time window, the 
fitting parameters are correlated: a change of the critical 
point position in the fit can be compensated, within cer- 
tain limits, by a change in the power-law exponent. For 
this reason, one hardly find data, which test the critical 
decay for more than a two decade time-window. It is 
also important to remark that testing MCT in systems 
which are described by a somewhat artificial kinetic rule 
is particularly interesting because it allows us to probe 
the degree of universality of MCT prediction beyond the 
original context (the actual Newtonian or Brownian liq- 
uid dynamics) and the closure approximations in which 
they were originally derived. 



A. Persistence 

To explore those universal features of relaxation which 
are relevant for a comparison with MCT we focus on the 
time evolution of persistence <fi(t), i.e. the probability 
that a spin has never flipped between times and t. It 
is a convenient and widely used characterization of the 
dynamics of facilitated spin systems. The long-time limit 
of 4>{t) plays the role of the nonergodicity parameter and 
represents the fraction of permanently frozen spins 

$ = lim <j)(t). (23) 
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When $ = the system is a liquid, while for $ > 
ergodicity is broken and the system is a glass. 
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FIG. 2: Persistence vs time for several values of q both in the 
discontinuous (q = 0, 0.2 and 0.4) and the continuous region 
(0.52) of the phase diagram, Fig. [TJ for r = 10"^ and several 
temperatures, either above and below T c (q). The solid black 
lines show the theoretical prediction, Eq. (|A8|) . for the plateau 
heights for T < T c (q) while the dotted ones signal the critical 
plateau at T = T c (q), q < 1/2. For temperatures above the 
transition, the long time value of the persistence is zero. 

We have been careful enough to disallow initial sys- 
tem configurations with a fraction of permanently frozen 
spins that would induce dynamical reducibility problems 
in our simulations (see the discussion in Ref. [39j for de- 
tails concerning these issues) . To this end we have chosen 
a ternary mixture with a rather negligibly small fraction 
of spins with / = 4, typically r = 10~ 3 , and not too large 
values of q (that is a fraction of spins with / = 3). In 
fact, we observed that the system becomes too overcon- 
strained and reducibility effects begin to appear above 
q ~ 0.55, at least for the system sizes, N = 10 5 — 5 x 10 5 , 
we have explored in our numerical simulations. At each 
time unity, called a Monte Carlo step (MCS), N attempts 
to flip an spin are performed. Simulation results for the 
persistence <f>(t) are shown in Fig. [5] for various values of 
q, both above and below the glass transition line T c (q). 
We see that above T c (q), every spin is able to flip so 
that the long-time limit of persistence is zero, confirming 
theoretical expectations and the absence of dynamical 
reducibility effects. Whereas for temperatures inside the 
glassy phase, T < T c (q), ergodicity is broken and the 
persistence attains a finite plateau which is in excellent 
agreement with the value of <£> analytically computed in 
Eq. (|A8[) , We also observe that the height of the crit- 
ical plateau $ c (shown as a dotted line in the figure), 
which is finite for q < 1/2, moves towards zero as q ap- 
proaches 1/2, and remains zero above it. In the previous 



section we have calculated how this $ c is approached 
when T — > T c (q) and we now compare the exact results 
with the numerical observations. 
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FIG. 3: Asymptotic fraction of blocked spins $ (plateau 
height) as the temperature approaches T c for several values 
of r and q = 0.52 in the continuous region of the phase dia- 
gram. Points are the result of simulations while the lines are 
from Eq. (|A8|) . The relevant temperature interval, in which 
$ ~ e~ , is not accessible in our numerical simulation when 
r 0. Inset: the same for q = 0.2 for which the transition is 
discontinuous (notice the different labels in the y-axis). 

In the main panel of Fig. [3] we show the behavior of 
$(T) for small values of r, as described before. The simu- 
lations are able to correctly detect the long term plateau 
as long as it is not too low, otherwise the time to attain it 
is exceedingly large and beyond our current computation 
capabilities. Indeed, the lowest point for which $ was 
measured in Fig. [3] corresponds to almost 10~ 2 , that is 
way before the beginning of the crossover to a linear be- 
havior. The linear in e long-term behavior manifests itself 
at a much smaller scale, not accessible in our simulation. 
On the other hand, in the region where the simulation 
was actually performed, the points are indistinguishable 
from the r = case and closer to the e 2 law. Only when 
3>(t) attains a value of the order of r, and after an exceed- 
ingly long crossover, the correct e behavior is observed. 
It should be noticed that larger values of r could diminish 
these strong preasymptotic effects, at the expenses of in- 
creasing dynamical reducibility effects. This effect could 
be at the origin of the discrepancy with the value of the 
measured exponent near the continuous liquid-glass tran- 
sition [fi9| . For the discontinuous liquid-glass transition 
instead we find that the square-root anomaly is very well 
verified. 



B. Relaxation times and critical decay exponents 

By measuring <fi(t) at various T and q, see Fig. [51 one 
can extract more quantitative information that will be 
useful to check important typical signatures of MOT. 
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When q < 1/2, on approaching T c (q) from the liquid 
side, the persistence <f>(t) presents the characteristic two- 
step relaxation form and we can reasonably estimate the 
characteristic times as follows: 



(f>(t)dt, 



<j)(t)dt, 



(24) 



where t c is the time at which persistence crosses the crit- 
ical plateau, <j>(t c ) = $ c . Notice, that according to these 
definitions, r coincides with 7> when $ c = (i.e. when 
the glass transition is continuous). 
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FIG. 4: Data collapse of persistence data for several values of 
q versus the rescaled time, t/r, where r is the characteristic 
time associated with the a relaxation. We obtain a good 
collapse for data such that <f>(t) < <3> c and, as $ c — > when 
q — 5- 1/2, the collapse region decreases as well, disappearing 
above the dynamical tricritical point. The collapsed region 
is well described by a KWW stretched exponential, shown as 
solid lines for q = and 0.2. The stretched exponent is 0.5 
for q = and 0.35 for q = 0.2. 

With these definitions of relaxation times, and after 
rescaling time by the structural relaxation time r, we can 
check the existence of universal scaling properties of per- 
sistence in the late a-regime. For times in the range over 
which <f)(t) < $ c , we find a nice data collapse of persis- 
tence function once the time is measured in units of r, see 
Fig. |4j The range over which this time-temperature su- 
perposition principle holds shrinks when q — > l/2~ (and 
so $ c — > 0), and consequently no simple scaling form 
of relaxation data at various temperature can be found 
in this limit. The scaling function is well fitted by a 
Kohlrausch- Williams- Watts stretched exponential func- 
tion with an exponent that decreases as q — > 1/2. 

To make a more stringent test of MCT we now exam- 
ine in detail the critical dynamics on the glass transition 
lines, for which MCT predicts distinctive patterns of uni- 
versal critical behavior. When the system is relaxing ex- 
actly at the critical temperature T c (q), MCT predicts a 
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FIG. 5: Equilibrium relaxation at criticality for discontinu- 
ous (q = 0.4, T c ~ 0.86) and continuous (q = 0.52, T c ~ 1.72) 
ergodic-nonergodic transition: solid lines in the main panels 
are power-law fits with exponents a ~ 0.23 and 0.31, respec- 
tively. Inset: /3-relaxation time vs temperature difference 
T — T c . The solid lines are power-law functions with expo- 
nents l/2a and 1/a, where the value of a is obtained by the 
fits in the main panel. 
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FIG. 6: Symmetric departure from the critical power law be- 
havior for curves at temperatures equidistant from T c , for 
q — 0.52. The closer the temperature gets to T c , the more 
symmetrical the curves become and the longer it takes to de- 
part from the critical power law. An analogous behavior is 
observed for q = 0.2 as well. 



9 



<3 0.2 - 
0.1 




1 ■ 1 ■ 1 ■ 1 ■ 1 

0.2 0.3 0.4 0.5 0.6 

Q 

FIG. 7: Exponent a obtained from the critical relaxation (cir- 
cles) and from the characteristic time as the critical line is ap- 
proached from above (squares). The solid lines are power-law 
fits that vanish as q — > 1/2, consistently with a logarithmic 
relaxation at that point. 
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FIG. 8: Estimate of the exponent b for q — 0.2 and r = 10~ 3 . 
For this value of q, T c ~ 0.58 and one obtains 7 ~ 2.81, 
a ~ 0.25 and, using Eq. Q, b ~ 0.61. This predicted value is 
in very good agreement with the simulation, as shown by the 
solid line whose slope is 0.61. 



power law decay: 

0(i,T c (g))-* c (g)~ t-W, (25) 

irrespective of whether $ c (<z) is finite or zero. The main 
panel of Fig. [5] shows that this power law behavior is very 
well obeyed, either above and below q = 1/2. Slightly 
above or below T c , the deviations from the power law are 
also predicted by MCT, Fig. [5] for temperatures equally 
distant from Tc, the departure from the critical power law 
is symmetric [l2j and starts roughly at the same time. 
Accordingly with MCT, the exponent a(q) is intimately 
related, in a way that depends on the nature of the glass 
transition, to the power-law exponent describing the di- 
vergence of the characteristic time 77, as the plateau at $ c 
is approached (the so-called /3-relaxation regime): when 
temperature gets closer to T c (q) one has Tp ~ e _1 / 2a ( | ?) 
for the discontinuous transition, and Tp ~ e - 1 / a (i) for 
the continuous one. These two distinct behaviors are 
tested in the inset of Fig. [5l Notice that the straight 
lines in the inset for describing the data of /3-relaxation 
time, measured with Eq. ((24]) . are not fitting functions 
but use the exponent a of the critical decay law measured 
in the main panel. The estimated r behaves as a power- 
law of e near T c and, more importantly, the exponents 
perfectly agree with MCT, for both the hybrid and con- 
tinuous glass transition. In Fig. [7J we collect the results 
of these two independent ways to measure the exponent 
a, for several values of q. They show excellent agreement. 
It is also interesting to observe that the exponent a de- 
creases when q — > 1/2, from both side of the continuous 
and discontinuous region. Although we were able to fit 
data of a{q) with a function that vanishes at q = 1/2, 
we were not able to approach this value of the exponent 
from the persistence at q = 1/2, at least with the sizes 
and times available to our computer simulation. 

Finally, we test the relation among the various expo- 



nents a, b and 7. For q = 0.2, the critical temperature is 
T c ~ 0.58, and we obtain 7 ~ 2.81 and a ~ 0.25 from the 
measure of r and t p , respectively, Eq. ©. Thus, Eq. ([7]) 
gives b — (27 — l/a) -1 ~ 0.61. Plugging this value of 
b in Eq. ([5J, we can compare the von Schweidler's law 
against the numerical data of simulations, Fig. [S] We 
observe that in the late /3 relaxation regime the agree- 
ment between the two is quite good. 



VI. DISCUSSION AND CONCLUSIONS 

We have extended the Fredrickson- Andersen facilita- 
tion approach to glassy dynamics by introducing spin 
mixtures with different types of facilitation rules [44|. 
Depending on the cooperative nature of the facilitation, 
they exhibit a crossover from a discontinuous to a con- 
tinuous glass transition. We have carefully tested by nu- 
merical simulation this approach on a Bethe lattice by 
analyzing the equilibrium dynamics near the liquid-glass 
transition lines. The results are fully consistent with pre- 
dictions of MCT and suggest that the cooperative facili- 
tation approach is the simplest microscopic on lattice re- 
alization of MCT. The existence of such a correspondence 
is by no means obvious because MCT is formulated for 
many-particle systems obeying Newtonian or Brownian 
dynamics, while the facilitation approach concerns par- 
ticles on lattice obeying an abstract-though physically 
motivated-kinetic rule. The close relation between MCT 
and systems with quenched disorder [2l|, [22J , suggests in 
fact that the dynamical facilitation description of glassy 
dynamics is not in conflict with the RFOT approach. The 
reason underlying this correspondence (that in a special 
case can be exactly proved |4(|) seems to be quite gen- 
eral [45j], even though it is not trivial to find the explicit 
expression of the mapping. To explore further the cor- 
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respondence between facilitated and disordered systems 
it would be certainly relevant to investigate the possibil- 
ity of inferring the exponent A, Eq. 0, from the MCT 
kernel or from the associated "free energy" as proposed 
in Refs. [10, Hl|- Also, it is interesting to observe that, 
while MCT struggles in the real space high dimensional 
(mean- field) limit [2!| [3(| , no such problem is observed 
in the Bethe lattice. 

The nature of the higher-order glass singularity and 
other peculiar dynamical features we observed is primar- 
ily controlled by the underlying bootstrap percolation 
process describing the formation of a spanning cluster 
of frozen spins. The geometric structure of this clus- 
ter can be fractal or compact depending on the nature 
of the cooperative dynamics. The connection with the 
bootstrap percolation problem can provide some clues 
on the fate of dynamical glass transition in finite spatial 
dimensions. Although we do not undertake this delicate 
task here, we mention that the existence of facilitated 
models with a dynamical glass transition in finite dimen- 
sions [6^ - l63 |. naturally suggests that the MCT scenario 
outlined above for higher-order singularity could be ob- 
served also in this circumstance. Also, we expect that 
the dynamical properties near the continuous and dis- 
continuous liquid-glass transitions are similar to those 
appearing in the F13 scenario (that differs from the F12 
kernel, Eq. (Q| , by the presence of a cubic term instead of 
a quadratic one) because the nature of such transitions 
is essentially the same in the two schematic models (i.e. 
a degenerate and generic A2 singularity). 

Although we focused on Bethe lattice with fixed con- 
nectivity and random distribution of facilitation, it is 
worth to remark that one can consider an equivalent vari- 
ant, in which the Bethe lattice is diluted and the facili- 
tation is uniform (44j. For example, results qualitatively 
similar to those reported above can be obtained when 
fi = 2 on every site and the local lattice connectivity Zi 
is distributed according to 

V{ Zi ) = (1-q) 5 ZiA + {q-r) 6 Zi , 3 + r 5 Zi}2 . (26) 

In this case, less connected sites turn out to have a 
smaller probability to flip, just as if they were less fa- 
cilitated. While the former variant has the advantage of 
being more easily implemented in finite dimensional sys- 
tems with fixed coordination number the latter is gener- 
ally more suitable for systems in which the local mobility 
does also depend on the geometrically disordered local 
environment. Thus, the crossover between the two glass 
transitions is obtained by varying either the local con- 
nectivity or the facilitation strength, and corresponds to 
a passage from bootstrap to standard percolation transi- 
tion. Whether the fractal cluster of frozen particles can 
be assimilated to a gel state is an interesting issue that 
it would be worth to address in a future work [H, [6(| . 
Finally, one of the advantages of the present approach is 
that one can quantify the contribution given by differ- 
ent types of facilitation to the glassy features by study- 
ing the dynamics of each species of spins separately. In 



this respect, it would be interesting to analyze and com- 
pare dynamical heterogeneities close to the continuous 
and discontinuous liquid-glass transition. 
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Appendix A: Calculations for determining the phase 
diagram and the fraction of frozen spins of 
facilitated spin mixtures 

We report here the detailed calculations to determine 
the liquid-glass transitions and the fraction of perma- 
nently frozen spins for models of facilitated ternary mix- 
tures. With the notation t + m + n we denote a ternary 
mixture composed by spins with facilitation / = I, m 
and n. For the sake of simplicity, we refer to a Bethe 
lattice with branching ratio k = z — 1 = 3, where z is the 
lattice coordination. We begin with the facilitated spin 
system we focused in the main text and then consider 
some other, qualitatively different, representative cases. 

1. Ternary mixture 2 + 3 + 4 

For the ternary mixture 2+3+4 with facilitation dis- 
tribution 

V{h) = (1 - q)5 ha + (q- r)6 fit3 + rS fiA , (Al) 

the probability that a spin is frozen 1 — B obeys the 
equation: 

1 - B = p [l - B 3 - 3(1 - q)(l - B)B 2 ] . (A2) 

This equation is always satisfied by B = 1, while an 
additional solution with B < 1 is obtained by solving the 
quadratic equation 

(3g - 2)B 2 +B + l-l/p = 0, (A3) 

which gives 

1 _ d _ 3(l-2g) + v /l-4(3g-2)(l-lM 

2(2 -3g) ' 1 ' 

The positive sign of the square-root is chosen because B 
is positive and should decrease as the temperature goes 
to zero. The continuous transition takes place when the 
above solution joins the B = 1 one. This happens when 

PM = i (A5) 
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Whereas the vanishing of the square-root argument de- 
termines the discontinuous transition: 



Pc(q) 



12q 



9 - 12q 



(A6) 



Using the relation between the temperature T and the 
probability p 



1 . P 

1 1 — p 



(A7) 



they provide us with a phase diagram showing a "tricrit- 
ical" line separating two glass transition surfaces. Let us 
focus first on the mixture 0+2+3, with 

T><Ji) = (1 - 9)Sf u o + (q- r)S fi>2 + rS fu3 . (A10) 

The probability 1 — B that a spin is frozen in this case 
satisfies the equation 

1 - B = p [q(l - B 3 ) -(q- r)(l - B)3B 2 ] . (All) 

Beyond the trivial solution B = 1 one also finds the extra 
solution 



one finally gets the dynamical glass transition lines T c (q) 
reported in the main text, Eq. (|2"T|). The crossover be- 
tween the continuous and discontinuous transition is ob- 
tained when the two expressions of T c (q) in Eq. ([2Tj) be- 
comes equal. This happens at q = 1/2 and corresponds 
to the tricritical temperature, Eq. (|22|) . 

The fraction of permanently frozen spins is 

$ = p[{l-B) 3 (A + 3B) + GqB 2 (l-B) 2 +Ar{l-B)B 3 ] 

(A8) 

+ (1 - p) [h 3 (4: - 3h) + 6qh 2 {l - h) 2 + Arh(l ~ h) 3 ] , 
where B is given by Eq. (|A2j) and h is 

h = p(l- B) [(1 - B){1 -B + 3qB) + 3rB 2 ] . (A9) 
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FIG. 9: Phase diagram, temperature T vs fraction of con- 
strained spins q, for the ternary mixture 0+2+3 (and 0+2+4) 
on a Bethe lattice with z — 4 and facilitation as in Eq. (|A10l l . 
The dotted/full lines are discontinuous/continuous liquid 
glass transitions for various r. The line passing through the 
circles is the tricritical line. 



2. Ternary mixtures 0+2+3 and 0+2+4 

Next we consider the diluted ternary mixtures of type 
0+2+3 and 0+2+4. These cases are interesting because 



1 - B 



6r - 3q + yjq 2 - 4(g - l/p)(3r - 2q) 



(A12) 



2(3r - 2q) 

As mentioned before, the continuous transition occurs 
when the above solution joins the B — 1 one, what hap- 
pens for p c = l/3r; while the vanishing of the square- 
root argument in the latter determines the discontinuous 
transition: 

4 2q - 3r 



Pc(q,r) 



3q 3q — 4r 



(A13) 



Using the relation between p and T one finally gets the 
dynamical glass transition surfaces, T c (g,r): 



T c (q,r) 



hi 



hi- 



4(2<? - 3r) 



12r(l - q) + q(9q - 



(discontinuous); 



1 



(continuous) . 

3r — 1 

(A14) 

Notice that the continuous transition exists in the inter- 
val 1/3 < r < 2/3 and does not depend on q. Interest- 
ingly, for any r in the range 1/3 < r < 1/2, there is a 
dynamical tricritical point whose location in the phase di- 
agram depends on the fraction of free spins. The location 
of this tricritical line separating the two glass transition 
surfaces is obtained by setting q = 2r, what gives 

77^TT=ln^-^. (A15) 
Ttric(g) 3g-2 

Sections of the phase diagram are illustrated in the Fig. [9] 
Since the equation term related to the spins with / 
. i is identical to that of spins with / = 4, the phase 
diagram of the mixtures 0+2+3 and 0+2+4 is just the 
same. While the related fractions of frozen spins differ: 
near the continuous transition line they vanish with a 
different critical exponent, $ ~ with (3 — 2 for the 
the mixture 0+2+3 and /3 = 1 for the 0+2+4 one. The 
phase diagram structure for mixtures 1+2+3 and 1+2+4 
is qualitatively similar. 

We report here only the fraction of frozen spins for the 
mixture 0+2+3, which is 

$ = pq(l- B) 3 (l + 3B) + 6prB 2 (l~ B) 2 (A16) 
+ q(l - p)h 3 (4 - 3/i) + 6r(l - p)h 2 (l - h) 2 , 

where B is given by Eq. (jAllI) and 

h = p(l - B) 2 [3rB + q(l - B)] . (A17) 
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FIG. 10: Phase diagram for the ternary mixtures 0+1+2 on 
a Bethe lattice with z = 4 and facilitation as in Eq. I|A18|) . 
The glassy region is located below the lines which represent 
discontinuous liquid-glass transitions for various r. 
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FIG. 11: Phase diagram for the ternary mixture 0+1+3 (and 
0+1+4) on a Bethe lattice with z — 4 and facilitation as in 
Eq. (|A20|I . The glassy region is located below the lines which 
represent continuous liquid-glass transitions for various r. 



3. Ternary mixtures 0+1+2, 0+1+3 and 0+1+4 

Finally, for completeness, we consider the ternary mix- 
tures of type 0+1+2, 0+1+3 and 0+1+4. The calcula- 
tion proceeds in a manner similar to those outlined above 
and, as expected, we find only one type of glass transition 
in these mixtures. We report here only the expression of 
their critical lines. For the 0+1+2 mixture with facilita- 
tion distributed as 

V(Ji) = (1 - «)*/«.o + (Q- r)S fi ,i + rS fu2 , (A18) 
the glass transition is discontinuous and we get 



1 



= - In 



T c (q,r) 

For the 0+1+3 mixture with 



r - 1 



(2q - 3r) 2 



4q - 12r 



(A19) 



V{fi) - (1 - q)S fi , + (g- r)6 fiA + rS fu3 , (A20) 



we obtain the continuous glass transition 



1 



T c {q,r) 



ln(4r 



!)• 



(A21) 



The above expression is valid for the case 0+1+4 as 
well. Sections of the phase diagrams are reported in the 
Figs. [TOl and [TT1 We notice that the glassy phases gener- 
ated in the discontinuous and continuous transition are 
both pretty stable against dilution (or random damage 
of the underlying lattice) . The fraction of free spins can 
be used in this type of systems as a way to control the 
plateau height. The phase diagram of binary mixtures 
can be easily obtained from the previous calculations by 
setting one of the weights to zero. 



[1] K. Binder and W. Kob, Glassy materials and disordered 
solids (World Scientific, Singapore, 2011). 

[2] W. Gotze and L. Sjogren, Rep. Prog. Phys. 55, 241 
(1992). 

[3] W. Gotze, Complex dynamics of glass-forming liquids 

(Oxford University Press, Oxford, 2009). 
[4] S. P. Das, Reviews of Modern Physics 76, 785 (2004). 
[5] D. R. Reichman and P. Charbonneau, J. Stat. Mech. p. 

P05013 (2005). 

[6] W. Kob and H. C. Andersen, Phys. Rev. Lett. 73, 1376 

(1994) . 

[7] W. Kob and H. C. Andersen, Phys. Rev. E 51, 4626 

(1995) . 

[8] K. Dawson, G. Foffi, M. Fuchs, W. Gotze, F. Sciortino, 
M. Sperl, P. Tartaglia, T. Voigtmann, and E. Zaccarelli, 



Phys. Rev. E 63, 011401 (2000). 
[9] K. N. Pham, A. M. Puertas, J. Bergenholtz, S. U. Egel- 
haaf, A. Moussai'd, P. N. Pusey, A. B. Schofield, M. E. 
Cates, M. Fuchs, and W. C. K. Poon, Science 296, 104 
(2002). 

[10] T. Eckert and E. Bartsch, Phys. Rev. Lett. 89, 125701 
(2002). 

[11] F. Sciortino and P. Tartaglia, Adv. Phys. 54, 471 (2005). 

[12] V. Krakoviack, Phys. Rev. E 75, 031503 (2007). 

[13] K. Kim, K. Miyazaki, and S. Saito, J. Phys.: Condens. 
Matter 23, 234123 (2011). 

[14] J. Kurzidim, D. Coslovich, and G. Kahl, J. Phys.: Con- 
dens. Matter 23, 234122 (2011). 

[15] E. Leutheusser, Phys. Rev. A 29, 2765 (1984). 

[16] U. Bengtzelius, W. Gotze, and A. Sjolander, J. Phys. C 



13 



17, 5915 (1984). [43 

T. R. Kirkpatrick, Phys. Rev. A 31, 939 (1985). 

A. Montanari and G. Semerjian, J. Stat. Phys. 124, 103 [44 

(2006). 

E. Siggia, Phys. Rev. A 32, 3135 (1985). [45 
S. P. Das, G. F. Mazenko, S. Ramaswamy, and J. Toner, [46 
Phys. Rev. A 32, 3139 (1985). 
T. R. Kirkpatrick and P. Wolynes, Phys. Rev. A 35, 3072 [47 
(1987). 

T. R. Kirkpatrick and D. Thirumalai, Phys. Rev. B 36, [48 
5388 (1987). [49 
T. R. Kirkpatrick, D. Thirumalai, and P. G. Wolynes, 
Phys. Rev. A 40, 1045 (1989). 
J. -P. Bouchaud, L. F. Cugliandolo, J. Kurchan, and [50 
M. Mezard, Physica A 226, 243 (1996). 

M. Mezard and G. Parisi, J. Chem. Phys. Ill, 1076 [51 
(1999). 

P. Mayer, K. Miyazaki, and D. R. Reichman, Phys. Rev. [52 
Lett. 97, 095702 (2006). 
S. M. Bhattacharyya, B. Bagchi, and P. G. Wolynes, [53 
Proc. Natl. Acad. Sci. 105, 16077 (2008). [54 
A. Andreanov, G. Biroli, and J.-P. Bouchaud, Europhys. 
Lett. 88, 16001 (6pp) (2009). [55 
R. Schilling and B. Schmid, Phys. Rev. Lett. 106, 049601 
(2011). [56 
A. Ikeda and K. Miyazaki, Phys. Rev. Lett. 106, 049602 
(2011). [57 
S. Franz, F. Ricci-Tersenghi, T. Rizzo, and G. Parisi, 
Eur. Phys. J. E 34, 102 (2011). [58 
L. F. Cugliandolo and J. Kurchan, Phys. Rev. Lett. 71, 
173 (1993). [59 
L. F. Cugliandolo and J. Kurchan, J. Phys. A : Math. 
Gen. 27, 5749 (1994). [60 
S. Franz and M. Mezard, Europhys. Lett. 26, 209 (1994). 
S. Franz and M. Mezard, Physica A 210, 48 (1994). 
L. F. Cugliandolo, J. Kurchan, and L. Peliti, Phys. Rev. [61 
E 55, 3998 (1997). [62 
L. F. Cugliandolo, J. Phys. A: Math. Theor. 44, 483001 
(2011). [63 
G. H. Fredrickson and H. C. Andersen, Phys. Rev. Lett. 
53, 1244 (1984). [64 

F. Ritort and P. Sollich, Adv. Phys. 52, 219 (2003). [65 
J. P. Garrahan and D. Chandler, Annu. Rev. Phys. 
Chem. 61, 191 (2010). [ 
J. Kurchan, L. Peliti, and M. Sellitto, Europhys. Lett. 
39, 365 (1997). 

M. Sellitto, J. Phys.: Condens. Matt. 14, 1455 (2002). 



M. Sellitto, G. Biroli, and C. Toninelli, Europhys. Lett. 
69, 496 (2005). 

M. Sellitto, D. de Martino, F. Caccioli, and J. J. Arenzon, 
Phys. Rev. Lett. 105, 265704 (2010). 
M. Sellitto (2012), arXiv: 1206. 2585. 

L. Foini, F. Krzakala, and F. Zamponi, J. Stat. Mech. p. 
P06013 (2012). 

J. Chalupa, P. L. Leath, and G. R. Reich, J. Phys. C: 
Solid State Phys. 12, L31 (1979). 
J. Adler and U. Lev, Braz. J. Phys. 33, 641 (2003). 
P. De Gregorio, A. Lawlor, and K. A. Dawson, in En- 
cyclopedia of Complexity and Systems Science, edited by 
R. A. Meyers (Springer New York, 2009), pp. 608-626. 
S. N. Dorogovtsev, A. V. Goltsev, and J. F. F. Mendes, 
Phys. Rev. Lett. 96, 040601 (2006). 

J. M. Schwarz, A. J. Liu, and L. Q. Chayes, Europhys. 
Lett. 73, 560 (2006). 

A. Montanari and G. Semerjian, Phys. Rev. Lett. 94, 
247201 (2005). 

N. S. Branco, J. Stat. Phys. 70, 1035 (1993). 

G. J. Baxter, S. N. Dorogovtsev, A. V. Goltsev, and 
J. F. F. Mendes, Phys. Rev. E 82, 011103 (2010). 

D. Cellai, A. Lawlor, K. A. Dawson, and J. P. Gleeson, 
Phys. Rev. Lett. 107, 175703 (2011). 

M. Sellitto and J. J. Arenzon, Phys. Rev. E 62, 7793 
(2000). 

Y. Levin, J. J. Arenzon, and M. Sellitto, Europhys. Lett. 
55, 767 (2001). 

H. C. M. Fernandes, J. J. Arenzon, Y. Levin, and M. Sel- 
litto, Physica A 327, 94 (2003). 

A. J. Moreno, S. H. Chong, W. Kob, and F. Sciortino, J. 
Chem. Phys. 123, 204505 (2005). 

F. Caltagirone, U. Ferrari, L. Leuzzi, G. Parisi, F. Ricci- 
Tersenghi, and T. Rizzo, Phys. Rev. Lett. 108, 085702 
(2012). 

G. Parisi and T. Rizzo (2012), arXiv:1205.3360. 

C. Toninelli, G. Biroli, and D. S. Fisher, Phys. Rev. Lett. 
96, 035702 (2006). 

M. Jeng and J. M. Schwarz, Phys. Rev. Lett. 98, 129601 
(2007). 

H. Ohta and S. Sasa (2011), arXiv:1112.0971. 

E. Del Gado, A. Fierro, L. de Arcangelis, and A. Coniglio, 
Phys. Rev. E 69, 051103 (2004). 

E. Zaccarelli, J. Phys.: Condens. Matter 19, 323101 
(2007). 



